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Abstract 

We study effective theories of an axion in spontaneously broken supersymmetric 
theories. We consider a system where the axion supermultiplet is directly coupled to a 
supersymmetry breaking sector whereas the standard model sector is communicated 
with those sectors through loops of messenger fields. The gaugino masses and 
the axion-gluon coupling necessary for solving the strong CP problem are both 
obtained by the same effective interaction. We discuss cosmological constraints on 
this framework. 



1 Introduction 



The hierarchy problem pQ and the strong CP problem [2] are two big misteries unexplained in 
the Standard Model of particle physics. Among various candidate solutions proposed so far, 
supersymmetry (SUSY) and the Peccei-Quinn (PQ) [3j mechanism are promissing ones. Both 
hypotheses rely on broken symmetries. The construction of the whole picture requires model 
building efforts of the symmetry breaking sectors and the mediation mechanisms to the Standard 
Model sector. 

The consistency needs to be checked when we combine SUSY and the PQ mechanism. If 
we are meant to solve the strong CP problem, it is not a good idea to assume a mechanism 
which causes another CP problem such as the SUSY CP problem. This gives a constraint on 
the SUSY breaking/mediation sector. Also, by SUSY the axion field is necessarily extended 
to an axion supermultiplet, which includes a CP-even axion field (saxion). The field value of 
the saxion (corresponding to the axion decay constant) needs to be stabilized, that requries a 
coupling to the SUSY breaking sector. The decay constant f a is constrained by cosmology and 
astrophysics. The viable region, 10 9 GeV < f a < 10 12 GeV, is called the axion window [3]. 

There are numbers of approaches for a supersymmetric axion model in field theories and 
string theories [5]-[2H]. The consistent model building tends to be rather complicated because 
one needs to carefully glue two symmetry breaking sectors. Especially, field theoretic approaches 
typically find new particles in both the PQ breaking sector and the SUSY breaking sector. The 
discussion of the viability is limited to specific models in such cases. On the other hand, the 
string theory approach tends to predict a too large decay constant [H] for the PQ symmetry 
breaking, such as the Planck scale, since that is the only scale in the theory unless there is 
some non-perturbative effects or a large compactification volume. The cosmological difficulty 
of the saxion (or moduli) field has also been pointed out |10 | 125 } [26] . The coherent oscillation 
of the saxion field and its decay produces dangerous particles like gravitino. There are severe 
constraints from the Big Bang nucleosynthesis (BBN) \27\ [28] and also the matter density of the 
universe |29j if the gravitino is the lightest supersymmetric particle (LSP). 

In order to avoid the complexity of the models for the discussion, we study simplified 
models of SUSY and PQ breaking in this paper. We construct effective theories which cover 
various microscopic models, although we do make a few assumptions motivated by phenomenol- 
ogy/cosmology and minimality. First, we assume the gauge mediation mechanism for the 
transmission of the SUSY breaking. This is a good starting point to avoid the SUSY FCNC 
and CP problem. We assume a direct coupling between the SUSY and the PQ breaking sectors 
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so that the saxion can obtain a large enough mass to avoid the cosmological difficulty. Finally, 
we assume that the same messenger fields (or the same effective operator) to communicate the 
SUSY breaking and PQ breaking to the Standard Model sector. 

In the next Section, we will propose a general framework for the simplified model. In 
Section 3, we will study properties of three representative models classfied by ways to break 
an i?-symmetry. The cosmological constraints are discussed in Section 4. In Section 5, we will 
give brief comments on the string theory approach. Section 6 will be devoted to the conclusion. 
In the appendices, we discuss supergravity corrections to the mass spectrum and list relevant 
computations of the decay rates of heavy fields. 



2 Low energy effective Lagrangian 



We construct a model where the PQ symmetry is non-linearly realized whereas SUSY is linearly 
realized. We introduce an axion chiral superfield A and a SUSY breaking chiral superfield X 
whose .F-component obtains a non-vanising vacuum expectation value (vev). In general, X can 
carry a PQ charge qx- The superfields X and A transform under the global t7(l)pQ symmetry: 

X^e iqxdpc >X, A^A + iOpQ. (1) 

Here #pq is a transformation parameter. Without a loss of generality, we take qx > 0. 



The effective theory for the PQ breaking sector is given by 



\{A + A^ + ^{A + A^ + ... 



(2) 



No superpotential cannot be written only with A. The parameter /q is the decay constant of the 
PQ symmetry breaking. We will take the coefficient C3 as a dimensionless parameter of O(l). 
The SUSY breaking sector has the Kahler potential: 

a{X^X) 2 {X^Xf 



X ] X 



+ ••■ 



(3) 



where Ao is the typical mass scale of the SUSY breaking sector and the parameter a can be 
chosen to be a = ±1 by an rescaling of the field. Xq is a prameter of order unity. The sign of 
a controls whether X obtains a vev [30]. The i?-symmetry is assumed here so that the SUSY 
breaking model described below will be justified. In the above Kahler potentials, we took A and 
X be dimensionless. Note that their values should be limited to \X\ < 1 and \A\ < 1 for the 
validity of the effective theor>|^] 

In order for X vev not to dominate the axion decay constant, it is required that /o/Aq > (X) if qx / 0. 
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The SUSY breaking is described by the Polonyi model [251 EH E2]. The superpotential 

is [331 El ES]: 

W SB =^ (W + A Xe-^ A ). (4) 

Here the X superfield has an i?-charge 2. The first term is a constant which is related to a 
fine-tuning of the cosmological constant in supergravity, Wq ~ Mp\/y/S. The gravitino mass is 
HqWq ~ m 3 / 2 Mp l . Here Mpi = 2.4 x 10 18 GeV is the reduced Planck mass. Note that A — > oo 
can be the SUSY vacuum; in that case a minimum with a finite (A) is a metastable SUSY 
breaking vacuum. 

Finally, we write down a coupling between the X and A fields which is necessary for 
stabilizaing the saxion field: 



Kax = A 



2 

c,0 



Cl (A + ^t) + | (j4 + A t)2^ x t X ( 5) 



o 



(6) 



Again, ci, ca, di and e\ are parameters of order unity. A negative value of C2 give a mass term 
for the saxion. The overall scale is set by 

A C)0 = min[A ,/o]. (7) 

The choice of this overall scale is motivated from the discussion below. 

For Ao <C fo, the above interaction terms are generated through, for example, a loop of a 
heavy field in the SUSY breaking sector whose mass carries the PQ charg^] 

AgxtAloglAoe^l 2 . (8) 

This corresponds to above interaction terms with A Cj o ~ Ao- 

For Ao 3> fo, on the other hand, the following terms can be generated after integrating out 
heavy modes in the SUSY breaking sector: 

S^SX^X, (9) 

where S = fae A . This case corresponds to A Cj o ~ fo- It is possible that the loops of fields in 
the PQ breaking sector generate 1//q suppressed terms which connect X and A. In that case, 



2 We can consider a O'Raifeartaigh model [36] in the UV scale such like W > = fJ 2 e~ qxA X + ^e qA (j>i<j)2 + 
\e q ' A X<i>\. 
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it is natural to assume that the \X\ term is also generated with the same suppression factor. 
Therefore, A Cj o ~ Ao ~ /o- In any case, one can summarize the overall scale as in Eq. Q. 

In summary the effective theory we consider is 

K = -Kaxion + -KSB + &AX , 

W = Wsb- (10) 
We take all the dimensionless parameters to be of 0{1). 

2.1 Unified origins of the axion-gluon coupling and the gaugino mass 

In this subsection, let us consider the origin of an axion coupling to gluons: 

d 2 9ATrW a W a , (11) 



where W a is the gluon superfield. This term is necessary to solve the strong CP problem. Once 
^-component of A acquires a vev, the same interaction induces the gluino mass. Since the 
term XTrW a W a is forbidden by the /^-symmetry, the above interaction term is going to be the 
leading contribution to the gluino mass. 

This unification of the axion coupling and gaugino masses can often be seen in string models 
(see [12] for the string theoretic QCD axion with intermediate scale decay constanlj^]). There, 
we have a coupling at tree level 

l - J d 2 6 + A) TrW a W a + c.c. (12) 

Here \jg\ = l/g 2 -i$/8-K 2 . Now Im(A) can be identified with (a linear combination of) the QCD 
axion and {g 2 ./2)F A is a gaugino mass, which can be comparable to or less than the gravitino 
mass. 

In gauge mediation models [391 HOI EH H2] , the above coupling can be obtained after inte- 
grating out messenger fields whose mass carries a PQ charg^J 

I^Mess = M e-^ A ^, (13) 

where q^q, = PQ(*&) + PQ(ty) and ^ and & are messenger fields. After integrating out the 
messenger fields, we have the axion coupling to gauge multiplets, 

\ J d 2 9 (Jp - q^N^^j TrW a W a + c.c. (14) 



3 See also topics related with LARGE volume scenario |37l 1381 121] . 

4 Even though we have an additional /^-breaking operator W = A' A.Xe~^ qx+q ^^ A '^^! , it is irrelevant to mass 



spectra, so far as the following conditions are satisfied: M > \' AXe~ qxA and M F A > X'AF x e~ 
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This is a hadronic axion model [33]. Here N^y is the number of messenger multiplets^J Im(j4) 
is the QCD axion and qyyNyq,(a/4Tr)F A is the gaugino mass, which can be larger than the 
gravitino mass. 



2.2 Standard Form 

We now define a particular basis to proceed the discussion. First, one can eliminate A in the 
superpotential by a redefinition of X by 

X -> e qxA X. (15) 

Next, we can also eliminate the (A + A^)X^X term by an appropriate shift of A. The vev (A) 
in this basis is vanishing up to small i?-breaking effects we discuss later. With new definitions 
of parameters, the superpotential and the Kahler potential are given by 

W = ^(Wo + AX), (16) 



K = -Ka x i on + i-CsB + KaX, 



(17) 



where 



K ■ — f 2 

lv axion — J 



\(A + A^ + ^{A + A^ + 



(18) 



Ksb = A 2 



x , x a(XtX) 2 (gtjQ 



18x1 



+ ••■ 



(19) 



Kax = K 



~-{A + tffX^X - d(A + A^) a{X \ X)2 



e(A + A^) 



(x^xy 

18x1 



H 



and 



(20) 



A c = min[A,/]. 



(21) 



The gauge kinetic term is 



We introduced a parameter k that depends on the messenger mechanism. 



(22) 



5 Once we require q^^N^y should be an integer, a value that q^^N^q, — 1 is desirable considering the domain 
wall problem [44] [45] 
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One needs to check if the shift we performed is within the allowed range in the effective 
theory. In terms of the original parameters, the shift is given by 

( A o) = -^-l(^r-c), (23) 

where 



2c 2 2 \q x 



'ci\ 2 1 2 / A x 1 



c 2 J q x c 2 \^c,o 
For \{A )\ < 1, we need q x = for A » /. 

3 /^-breaking Models 



C=l/(-) + • (24) 



The scalar potential can simply be calculated from (16) and (17): 



V = K xtx \d x W\ 2 . (25) 

Recall that a small i?-symmetry breaking effect is required to produce gaugino mass; a vev of 
F A is vanishing without R- breaking effects. 

We consider three types of i?-symmetry breaking in the following. The first one (Model 
1) is a model with spontaneous i?-symmetry breaking, that can be achieved by taking the a 
parameter to be —1. The second model (Model 2) is to add a small explicit /^-breaking term in 
the Kahler potential: 

AK = e K A 2 {X +Jft)(A + At). (26) 

Finally, the third model (Model 3) is to add a small explicit /^-breaking term in the superpo- 
tentialB 

X 2 

AW = efi 2 A- — . (27) 

We assume a = 1 in Model 2 and 3. 

We discuss in the following the mass spectrum in each model. 



6 AK can be, for example, generated from an additional i?-breaking messenger term W = 
A'AXe~ (9x+9 *4 )A ** [46]. On the other hand, AW can originate from a non-renormalizable superpotential 
W = ^s 2 X 2 , where s = fe~ qxA . 
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3.1 Model 



Before going to the discussion of the i?-breaking models listed above, we comment on the effect 
of gravity mediation. Since the i?-symmetry is explicitly broken in the supergravity Lagrangian, 
F A is induced and the gauginos obtain masses of 0(m 3 / 2 ). 

The axino (the fermionic component of A) and the saxion (Ar) can also obtain masses. The 
axino mass is of 0(7713/2) arid 

m A R = X V6c J • ( 28 ) 

The mass of the scalar component of X is 

mx = V3 3/2 A . (29) 

This is the minimal model of the supersymmetric axion. However, the gravity mediation 
potentially has problems of large CP violation which we are trying to solve. Therefore, we 
consider cases with m 3/ / 2 <C 100 GeV where the supergravity effects is subdominant. 

3.2 Model 1: spontaneous i?-breaking model with a = — 1, AK = 0, AW = 

For a = — 1, X = is unstable and the R symmetry is spontaneously broken by its vev. The 
potential is stabilized through the Kahler potential K ~ — (X^X) 3 . The vevs of X and A are 

4c V A , 

(A) = 2 ~^xl (30) 

At the vacuum, the masses for the CP-even scalar fields are obtained to be 

m 2 - 2// m 2 _ 2 ^VM 2 r311 
Xr ~ A 27 ' Ar ~ \T J ' ^ ' 

Here X = (xq + X^)e* X/ , A = Ar + iAj and the above expressions are correct up to of O(xq). 
The parameter xo should satisfy a condition^] 

^0 < 1, (32) 

for the effective theory to be valid. Since the i?-symmetry is spontaneously broken, there is a 
nearly massless .R-axion Xj. The Xj field and the goldstino (the fermionic component of X) 



7 {X) ~ xq < 1 also satisfies the condition to obtain the stable vacuum. This is similar to Model 3. 



s 



acquire masses through a supergravity correction: 

2 2/i 4 /- 2 / MpA 

up to corrections suppressed by the Planck scale and of 0(xq). The axion Aj remains massless 
at this stage. 

The axino a, which is a fermionic part of A, obtains a mass via a Kahler term: 

K = -K 2 c -^{A + A^fX^X. (34) 



It is given by 



mi = ^ j (3S) 



and a fermionic part of X, ipx, is the goldstino which is absorbed into gravitino. Finally, the F 
component vev for the axion is 

?A _ (3d-2eVA^V 4 (36) 



6 p \A 
This gives the gaugino masses, My 2 = (ka/A / K)F A . 

3.3 Model 2: explicit i?-breaking model with a = 1, AK ^ 0, AW = 



to 0(ej v 



This model is i?-symmetric if we neglect ex- Thus we have the following scalar mass spectra up 

m L = m 2 Xl = jj, (37) 



-L = ^ ¥ ■ (38) 



o _ 2c> 4 /A, 

* " T 2 " V A 

Here we denoted as X = Xr + iXj. Of course, the axion Aj and goldstino ipx remain massless 
at this stage. 

The shift of the vevs at the leading order in the e# expansion ar^] 

8X = e 3 K C 3 ^, (39) 

M = ^^rA¥) 2 - (40) 



'2c/ 2 V A 



For small etc, the SUGRA correction to 8X may have the same order of the magnitude as the following 
result. 
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The mass mixing between Xr and Ar is generated by ex- 

ro ^ = T7 r ^" 6c UJ + a*J" ( } 

In order not for the vev shifts to be larger than 0(1) and also not to destabilize the vacuum by 
the mass mixing, we obtain a condition 

«<(£)'. m 



The F-component vev of A is: 

FA = ^. (43) 
This induces the gaugino masses and the axino mass through a Kahler term (f 2 C 3 /3l)(A + A^) 3 : 

m a = C 3 F A . (44) 

3.4 Model 3: Explicit .R-breaking model with a = 1, = 0, AW ^ 

In this model, scalar and goldstino mass spectra are the same as Model 2 up to 0(e 2 ). Here the 
vev shift of X is 

SX = -e, 5 A = 0. (45) 

Therefore, we have a condition to make our effective theory valid: 

e < 1, (46) 

for validity of the effective theory. Then the axino obtains a mass via K = A 2 c(A + A^) 2 X^X 
as in Model 1: 

m~ a = -ce.^A . (47) 



The F-component of A is 



f* - ^ (£ ) • («) 



10 



3.5 Summary of mass spectra and constraints on parameters in three models 



Here we summarize the mass spectra of three models in Table [T] Recall we have several 
parameters: 

c, d, e < 0(1), C 3 = O(l). 

There are also constraints on i?-breaking parameters in each model: 

Model 1 : x < 1, 
Model 2 : e K < 



A 



Model 3 : e < 1. 



(49) 
(50) 
(51) 





Model 1 


Model 2 


Model 3 


m x R 


(V2F X ) 2 




(F*) 2 


m 2 Xi 








2 

Ar 


^^2~cF x ) 2 


^^2dF x Y 


{^fV2cF x ) 2 


ma 


* c(^)V 


C 3 F A 




F A 


(2e-3d)x 3 hi vX 
6 P r 


-e K • ^F x 


2 



Table 1: A table of mass spectra in three models for a gauge mediation with light gravitino. 



Here AF X 



given by M 1/2 ~ k§F A and m ~ ^F^ 4 . 



-\/3m3/2Afpi and A c = min[A, /]. The soft mass in the visible sector are 



The masses in the table should be smaller than the cut-off scales, A/4ir in the SUSY breaking 
sector or / in the PQ breaking sector. Therefore, we have a constraint: 



A > V^7TjU. 

4 Cosmological constraints 



(52) 



In this section, we discuss cosmological constraints on three models defined above. We consider 
decays of the coherent oscillations of the saxion and Xr (the Polonyi field) and discuss constraints 
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from the BBN and the matter energy density of the universe. We assume that the LSP is the 
gravitino and the next-to-lightest supersymmetric particle (NLSP) is the bino. The bino NLSP 
gives more stringent constraint than a case for stau NLSP. 

In the following discussion, we fix the PQ breaking scale / to be 10 10 GeV, bino mass to 
be 100 GeV and all dimensionless parameters (except for xq, €k and e) to be of 0(1). Once we fix 
/ and the gaugino mass, the parameters we have are A and the gravitino mass (or equivalently 
the i?-breaking parameters in each model (xq, ex, or e)) aside from the 0(1) parameters. See 
Figure [l}j3] for the parameter regions constrained by the validity of the effective theory. 

In general, light scalar fields such as Xr and Ar are problematic for cosmology since their 
coherent oscillations and late-time decays would produce too large entropy and also produce 
unwanted particles. In the models we are discussing, the Polonyi and the saxion can be much 
heavier than the Standard Model particles due to the direct coupling to the SUSY breaking 
sector. This situation helps to cure the problem. We consider here the constraints on the model 
parameters from the successfull BBN and the matter density of the universe [10, 25, 26J. 

We define several important temperatures for discussion. The decay temperature T? is the 
one at which a scalar condensate cj>(= Xr,Xj,Ar) decays in the radiation-dominated universe. 
The domination temperature T^ ora is the one at which <fi dominates over the energy density of 
the universe [1711351115] : 



Here g*(T$) is the effective number of light particles at TJ, is total decay width of (p. Tr 



starts to oscillate and Acf) is the initial amplitude. These will be important for our discussion 
below since <$> does not dominate over the universe if T% > T^ om . 

4.1 Model 1 

We first calculate the decay widths of particles in Model 1. 
• Xr (Polonyi) and Ar (saxion) decay 




(53) 



(54) 



is a reheating temperature by an inflaton decay, T c ' 



-<t> ~ 

osc — 



O.3y/m0Mpi is a temperature when (j> 
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Through the interactions K = K 2 [X^ X — (X^ X) 2 /4], the Polonyi field decays into i?-axions or 
gravitinos. The total decay width is given by 



T Xr ~ F(X R -> XjXr) + T(X R -> ^3/2^3/2) (55) 

1 my 1 W-v 

+ 7^ 2"%^- (56) 



32ttx 2 A 2 96tt mj^M 2 ^ 
The vev xo can be expressed in terms of the i?-axion mass as 

x A = 2V3M Pl ( . (57) 

For the saxion, the main decay modes are Ar — > aa and Ar — > ^3/2 + 5, which originate 
from K = f 2 (C 3 /3)(A + At) 3 - A 2 (c/2)(A + A*)*X*X. The decay width is 



F Ar ~ r(^ -> aa) + r(,4 R -> V3/2 + 
_ C 3 2 , 1 m\ R 



a 



~ 32vr f 2 + 96n m 2 3/2 M 2 i (58) 

The axino (a) produced by the saxion decay subsequently decays into gravitinos or the Standard 
Model particles with the decay width: 

Ta ~ T(a^^ 3/2 + a) + T(a^ X + g) (59) 

1 m~ , T or m~ , , 

+ ^^3^f> (60) 



96vr m 2 /2 M 2 L 9 256^ 3 / 2 ' 

where iV fl (=12) is the number of the decay mode 

Just by looking at the main decay modes, one can see that it is problematic if the coherent 
oscillations of the Polonyi field or the saxion field dominate over the energy density of the 
The decay products are stable (or long-lived) particles such as gravitinos and axions 



universe 
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which contributes to the matter energy density of the universe. The overproduction of those 
particles needs to be avoided for viable cosmology. 

For of 0(Mp\) initial amplitudes of the Polonyi and the saxion, the condition for the matter 
density f2 ma tter ~ 0.2 requires a low enough reheating temperature after inflation such as Tr < 1 
MeV. Such a low reheating temperature may be barely consistent with the BBN. However, 
we do not consider this case since 0(Mp\) field values are beyond the validity of the effective 



9 There is also a decay mode a — > 1P3/2 + Xi, which can be the main decay mode in a narrow region in the 
parameter space. 

10 If the saxion s — fe is captured at the origin s = during the inflation, the thermal effect via messenger 
fields becomes relevant to the saxion potential. For such a case, the saxion can dominate the energy density of 
the universe fT71. 
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theory. One can consider possibilities that the initial amplitudes are at AX/j = 0(A) [51] or 
AAr = O(f) |20j, since A and / are the unique (cut-off) scales for X and Ar, respective!^] For 
such a case, the Polonyi (the saxion) decays before dominating the universe, with the parameter 
region that we are considering, since conditions T^ m < and T d R < T^ R are always satisfied. 
Hereafter we consider such a case to search for viable parameter regions. 

• Xj (i?-axion) decays 

The i?-axion decays into two gauginos if it is kinematically allowed. The interaction term in the 
Lagrangian is [19] 

Xj Mi to 

-i— — AA + c.c (61) 

^2x A 2 v ; 

Here A is the MSSM gauginos. The total decay width is 

T Xl ~ r(Xj^AA) (62) 
N g (M l/2 \ 2 ( 4M 1 2 / ^ 12 



For mxj < 2M 1 / 2 ) the channel Xj — > AA is closed, then Xj — > bb is the main decay mode 
through the mixing between Xj and the CP-odd Higgs boson (A) in the MSSA-Q The mixing 
is obtained through the Bfj,-term. Although we need a concrete model to generate i?^-term 
to discuss that interaction, there is always a contribution from a one-loop diagram with the 
gaugino mass insertion |53l 154] . even if we have B\i = at the messenger scale [55]. If that is 
the donimant contribution, there is an interaction term to mix Xj and A: 

i m 2 A sin 2/3 ... 

W2-^r Lx,H " Hd+c - c -' <64) 

where we have used a relation B[i = m 2 A sin 2/3/2. From this interaction and Yukawa coupling, 
the decay width is found to be |49| : 

r wv ,tn 3m X/ / m\ sin 2 m b V /" 4mf 

r^r(x,-wo * — XoA m 2_ m 2,J f-^- ^ 

For further smaller rax f , X/ — )• rf can be the main decay mode with large tan f3; should be 
replaced with m T . In the parameter region that we are interested, T, Xl > 1 MeVR 



A possiblity A.Ar = O(vToBgut) was also considered in Ref. |52) . 

12 We neglected a decay mode X/ — > tt via the similar interaction, since we assumed mt > nig = 100 GeV. On 
the other hand, note also that so long as tan ft > \Jmtjrrn, ~ 6.4 this mode is suppressed, compared to Xi — > bb. 
In a whole computation, we assumed vtia and mxj do not degenerate. 

13 If the decay of Xi is too late, it can influence the BBN. This excludes a small parameter region around 
A ~ 10 13 GeV and m 3/2 ~ 40 MeV in Figure 1 [27]. 
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On the other hand, because xoA is a normalization of the i?-axion, its initial amplitude is at 
most on the order of xqA. The temperature at which oscillating i?-axion dominates the universe 
is given by 

Here T*£ ~ 0.3yW^Wpi - 1-6 x 10 10 GeV (m X/ /1.2TeV) 1/2 is a temperature where the R- 
axion starts to oscillate. The relation m^^/mxj ^5x 1CP 4 holds for 7713/2 < 40 MeV which we 
need from the BBN constraint we discuss later. Then we obtain T^ 1 ^> T^JL in our model; the 
i?-axion does not dominate over the universe. 



• Bino NLSP and BBN constraint 



We have seen that the Polonyi field, the saxion, or i?-axion do not dominate over the energy 
density of the universe provided that the initial amplitudes of the Polonyi field and the saxion 
are of O(A) and of O(f), respectively. Then we need to consider the bino abundance since it is 
a long-lived NLSP, which can disturb the BBN. The thermal abundance of binos is given by [56] 



below the freeze-out temperature T? ~ m^/30. 

Since the binos decay into gravitino at a later time, there is constraint from the BBN on the 



lifetime as mentioned above. With the yield of Eq. (66), this can be translated to the bound on 
the gravitino mass which is |28j 

m 3/2 <40MeV for m § = 100 GeV. (67) 

For a larger bino mass or for the stau NLSP case, the constraint is relaxed, such as 7713/2 < 1 
GeV. See Figure [T] for the allowed region. 

There are contributions to the bino density from the decays of the heavy field such as Xr, 
Ar and a. Since they are heavy enough, the binos produced by those decays are thermalized 



and thus such contributions are already taken into account in Eq. ( 66 ) . The decays of Xj can 
be later than the freezing-out of the bino-pair annihilation. This non-thermal contribution is 
smaller than the thermal piece when 

Tr < 10 6 GeV, (68) 
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which is required later by the constraint from the gravitino abundance. 
In summary in Model 1 the gravitino mass is constrained to be 



1(T 4 GeV < m 3/2 < 40 MeV, 
and the cut-off scale (or dynamical scale) of the SUSY breaking sector is 

10 8 GeV < A < 10 13 GeV, 

for / = 10 10 GeV. 



(69) 



(70) 




L*S1Q[ ] 

lGeV 



Figure 1: A figure of allowed parameter region for M\/2 = 100 GeV. Recall that we have 
constraint that should satisfy ■s/4n/j, ^ L, xq < 1 and m 3 / 2 ^ 40 MeV. Here L = {\/2xq) 1 I 2 A. < A. 



Candidates of cold dark matter 



The thermal contribution to the gravitino density ^3/2 is |57j : 



Q th 

Ji 3/2 



9 T R \ /30MeV\ / M- g - 2 

' V 3.1 x 10 6 GeVy V m 3/2 J V600GeV 



(71) 
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Here Mg is gluino mass at TeV scale. Once one fixes gluino mass, the gravitino by thermal 
scattering can be dark matter of the universe for Tr < O(10 6 ) GeV x (to 3 / 2 /O(10) MeV). 

Note that we also have non-thermal component to the gravitino density. The important 
contribution is from the Polonyi decay: 

Y W * 2m^ R B ^- (?2) 
Here B 3 ,^ is a branching ratio of a decay mode Xr — > 2^ 3 / 2 - Then we obtain 

^-»- 23 (lo4v) S (io^v)^- < 73 > 

Here = O(l). For A ~ 10 13 GeV and m 3 / 2 — 40 MeV (near the upper right boundary of 
the Figure [TJ, the non-thermal component can be as important as the thermal one. Considering 
the BBN constraints, the contribution to the gravitino density from the decays of binos, saxions, 
is much smaller than this contribution for Tr < O(10 6 ) GeV and 7743/2 ;$ 40 



and the axinoi 



11 



MeV. This is because of low Td om and small branching ratios into gravitinos. 

The axion is also a candidate for dark matter [35]. The abundance is given by 



. \ 2 / f \ 7/6 

'-'mis \ / J a 

IT) V 10 12 GeV 



n.^i.4(^) It^) , (74) 



where O m i s is misalignment angle of the axion. 
4.2 Model 2 and Model 3 

Essentially, the discussion is parallel to Model 1. A difference is that the i?-axion now has a 
similar mass to the Polonyi field, and thus we do not need to consider it separately As in Model 
1, the domination of the Polonyi and the saxion fields would produce too much gravitinos, and 



14 Recently, a paper |68| discussed the constraint on the reheating temperature obtained from the matter 
energy density of the thermally produced axinos [BS] or gravitinos from their decays. In our model, the gravitino 
abundance from the thermally produced axino is given by 



~ Tr \( rn3 ' 2 \ ( w GeV V ^ 

" 3/2 -°- 2 v 8.7xl05GevH30MevM f a ) ^ 



Here we have replaced Mg with v / 6a3m3/2Mpi/(47r/ a ) and used (73 ~ 1 in eq. |7l| , and By 2 is a branching ratio 
of axino decay to gravitino. As one always finds B3/2/ fa ^ 10 -24 GeV -2 for f a < 10 12 GeV, this can give more 
stringent constraint on Tr together with eq.( |71[ ). The allowed region for the reheating temperature is 



T R <m&, or m a < Tr < 8.1 x 10"GeV 



f a \ /600GeV\ 



lO^GeVy V M- 9 ) ^ 2 
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therefore the initial amplitude should be small enough. With O(A) and 0(f) for the sizes of the 
initial amplitudes, the problem can be avoided. 

The viable parameter regions are shown in Figures [2] and [3] which are very similar to Figure [TJ 
Therefore, the viable range we obtain is the same as Eqs. ( |69[ ) and ( 70 ). We summarize in Tabled 
the numerical values of the masses and decay widths in the parameter range of our interest. 




l*siq[ ] 

lGeV 



Figure 2: A figure for an allowed parameter region for Model 2. This is similar to the Model 1 



except for ejc ^ V . Here we defined r\ 



and Ek = ckV 



5 Brief comments on axions in IIB orientifold/F-theory GUTs 

In string therory, besides the field theoretic axions, we often obtain very light string theoretic 
axions via moduli stabilization in the low scale SUGRA |58} 159] . In general, the number of 
axions is estimated as 

(The number of axions) = (The number of moduli fields) +1— (The number of terms in the W). 

Here W is the superpotential and a factor unity comes from the i?-symmetry. This is because 
PQ symmetries of (moduli) fields and the i?-symmetry produce candidates of the axion whereas 
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™3p. 
lGeV 

Figure 3: The similar figure for Model 3. This is similar to the Model 1 and 2 except for e < 1. 
The discontinuity of e originates from the fact F A oc d = (d\ + 2qx) and we took qx = 1 for 
/ > A while qx = for / < A. 
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Model 1 


Model 2 


Model 3 


m XR (GeV) 


10 5 - 10 s 


10 4 - 10 s 


10 4 - 10 8 


T Xr * T Xr ^ 3/2 (+r$£_2X,) (GeV) 


IQ-io _ 10 4 


10~ 16 - 10 4 


10" 16 - 10 4 


m Xl (GeV) 


5 - 5 x 10 3 


10 4 - 10 s 


10 4 - 10 s 


r v ~ rMl nr T M1 nr F M2 ' 3 fOeV) 


10-24 _ lo-n 


10~ 16 - 10 4 


10" 16 - 10 4 


mi fi (GeV) 


10 5 - 10 7 


10 4 - 10 s 


10 4 - 10 8 


Ta r ^ r AH _».2a + r Ajl ->V3/25 (GeV) 


10~ 7 - 10 


10" 10 - 10 


10~ 10 - 10 


m 5 (GeV) 


5 x 10 4 - 10 6 


10 4 


10 4 - 10 6 


Ta — Ta^ 3/2 a + ^a^Xg (GeV) 


10~ 12 - 10~ 5 


10-12 _ 10-10 


10~ 12 - 10~ 6 


ii-breaking parameter: (xq, e) 


10-0-8 _ i 


10~ 3 - 1 


10" 1 - 1 



Table 2: A table of examples of numerical values in each model. We took all dimensionless 
parameters as of order unity, = 100 GeV, / = 10 10 GeV and m 3 /2 ^5 40 MeV. Here we 

defined rj = A c /A. 



independent terms in the superpotential kill them, supposing the Kahler potential K preserves 
these symmetries. When this counting gets negative or zero, we do not have any light axions. If 
there are very small subleading terms violating PQ symmetries in W or K, they give very light 
mass to the axions. If we are to identify one of the axions with the QCD axion, the quality of the 
PQ symmetry needs to be checked for solving the strong CP problem: 5m 2 < 10 _11 (ma CD ) 2 . 
Here axion mass 5mi is a contribution from non- QCD effects and (m$ CD ) 2 is the QCD axion 
mass just from the instanton. 

For a string theoretic (QCD) axion, we often encounter a logarithm type Kahler potential like 
-Kaxkm = — Mp^ log(A+A), while we can see also quadratic Kahler potential -ftTaxion = ^f(A+A*) 2 
near the singularity. For logarithm case, by the expansion around the vev like the standard form, 
we can obtain the quadratic one with a decay constant / = ^^ A t)) ~ /* = 0(M s trmg), where 
-^string is the string scale. Then we can also obtain K = (X* X) 2 / f 2 after integrating out a 
massive gauge boson of (non-) anomalous f7(l)pqj^j Its contribution to the X mass is on the 

For non- anomalous U(l) the gauge boson mass is lower than the string scale bacause a matter-like field always 
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order of F / f. This corresponds to the case with A ~ /, once we could manage to obtain 
/ at an intermediate scale. One should also consider the cosmological constraints since there 
typically exist light moduli in a low scale SUSY breaking scenario [38J. (See also a recent paper 

ED].) 

In Ref. [33], it has been studied a model which has the same superpotential as in Eq. Q. 
There, C/(1)pq is an anomalous gauge symmetry at a high energy scale, and X is the Polonyi 
field. A is a string theoretic axion multiplet and is describing the 4-cycle volume on which a D7- 
brane holding the C/(1)pq is wrapping. In the paper (X) = f a /V2 is obtained and A is absorbed 
into the U(1)pq gauge multiplet via D-term moduli stabilization since we have f a <C /a- Here 
/a is a decay constant of A. Hence, at low energy, the relevant field is only the Polonyi field, 
X; the field Arg(X), the i?-axion, is identified as the QCD axion. However, from the discussion 
in Refs. |2U| and |61j . the axion becomes either massive or the decay constant is unacceptably 
large. Therefore we need to modify the superpotential to have the QCD axion successfully. 

For this purpose, we may consider an effective superpotential on intersecting seven branes 
induced by stringy instantons |62| . for instance: 

W = eM Pl sXe~ qBB + X'(s- s' -e'M^ ie - qBB ) + s'^ + ■■■ (75) 

with ?7(1)pq x [/(I)pq' symmetry where both U(l) symmetries are anomalous. Here and ^> 
are messenger fields, dots represent constant Wo, heavy fields in the SUSY breaking sector or in 
the PQ breaking sector, the moduli stabilization sector and so on. X' is a Lagrange multiplier 
superfield. e and e' are complex structure moduli/dilaton and other moduli contributions 16 
besides B. B is now a string theoretic axion multiplet describing a proper 4-cycle which would 
not intersect with a D7-brane holding [7(1 )pq (or its orientifold image brane), and transforms 
non-linearly under £7(1) PQ /; we may find PQ(e- qBB ) = PQ(X') = 0, PQ(X) = q x = -q^ 
and q B = PQ'(s) + PQ'(X) = PQ'(X') ^ 0. s and sf may appear as s = fAe~ qxA and 
s' = fA£ qxA through a constraint where we defined Mp\(e'e~ qBB ) = f\ = 0(f%) and we obtain 
X' = via an F-term equation. Stabilization of these fields, however, can depend on the model. 
For instance, let £[/(i) PQ and Cu(i) PQ , be Fayet-Iliopoulos (FI) terms in [7(1)pq and [7(1)pq/ 
respectively. Note S,u(i) PQ , includes B while Cu(i) PQ can depend on other modulus, say C. Then 
D-term potential (with non-trivial world volume flux) would almost fix C and B at Cu(i) P q — 

cancels its Fayet-Iliopoulos term; in general we obtain K = (X" 1 X) 2 /m^j^ , where is the (non-)anomalous 

U(l) gauge boson mass. 

16 We will assume e and e' are singlet under two U(l)s for simplicity below, and they might become important 
to get proper scales. On the other hand, instead of e'e _9BS , we could have an another moduli contribution 
e " e - q B' B ' (£' B) for instance. 
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and Cc(i)p Q / — 0) leaving X and A light mode^J Thus each U(l) gauge multiplet would get 
massive by eating B and C, assuming decay constants of B and C are of 0(M s t r ing) which are 
larger than /a and |(-X")|. Hence we would have the global l7(1)pq symetries in the low energy 
scale 18 , though another global C/(1)pq/ will be broken by a stringy instanton. Then we could 
gam W = W + {x 2 Xe~ qxA in the low energy. Further analyses would be beyond our scope since 
they depend on string model including issues of string moduli stabilization (see for instance 

m\). 



Finally, we would like to mention moduli stabilization in type IIB orientifold models. For 

moduli mediation scenario, we can get (more than) TeV scale moduli masses. This is because 

moduli masses are typically related with SUSY breaking gravitino mass because of axionic 

shift symmetries. They should be heavy enough to avoid cosmological problems. For gauge 

mediation scenario, we may need light gravitino mass and may get subsequently light moduli. 

Hence, for such a case with a field theoretic QCD axion, supersymmetric moduli stabilization 

which breaks the shift symmetry should be considered, such that we obtain heavy moduli: 

h + 

(W) » 0, (d moduh W) w via fine-tuning [63]. For instance, W = W +J2i M (Aie _aiTi +Bie _&iTi ) 
with a fine-tuning of fluxes Wq < 1 or non-geometric flux compactifications [63] might be viable. 
To obtain large volume, T» S> 1, we will need further fine-tunings among Aj, Bi,a,i, and hi. 
The remaining Kahler moduli may be stabilized via D-terms as mentioned above. Then we 
get _p modul1 = and much smaller constant term in the superpotential to obtain light gravitino 
mass, and moduli mass are decoupled to the gravitino mass: (W) <C (<9^ oduli W0 in the Planck 
unit. 



6 Discussion and conclusion 



We discussed supersymmetric effective theories of the axion field and the Polonyi field, which 
can be charged under global J7(1)pq symmetry. We assumed that the SUSY breaking and the 
PQ breaking sectors are directly coupled whereas the visible sector is communicate with those 
sectors through messenger fields. For a concreteness, we construct theories by Taylor expansion 
in X and A. To compute those in a simpler way or as generalizations, we transformed theories 
to "the standard form" with (partial) C/(1)pq transformation after expansion around the saxion 

17 We assumed an instanton depending on C does not have two fermionic zero modes related with Af = 1 
supercharge broken by the instanton or its coefficient in the superpotential has a vanishing vev; dcK — could 
mean the potential minimum and 8bKW + M^ObW = would be simultaneously satisfied there, though £ < f\ 
would be sufficient for us to get a proper axion decay constant. 

18 We could have another possibilities without using an anomalous t/(l) gauge symmetry. For accidental axion 
models in heterotic orbifold, see [18]. For a discrete i?-symmetry argument which heterotic string models might 
have, see [23] . 
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vev {Aq). Then we considered three ii-breaking models to obtain large gaugino mass in a gauge 
mediation with light gravitino mass such like 7713/2 ^ 40 MeV to avoid BBN constraints. In 
those models, axion multiplet plays important roles in mediating axion and SUSY breaking to 
the visible sector. Both the Polonyi and saxion can be heavy because the former has the low 
scale cut-off A and the latter has a direct coupling to the Polonyi field. However, they must 
not dominate the universe after the inflation. If the initial amplitude is of the order of the 
cut-off scale of the effective theory, such dominations can be avoided. For Tr < O(10 6 ) GeV or 
Tr < ma, the axion or gravitino can be dark matter of the universe. 

We have several things we did not consider here explicitly. For instance, since inflaton can 
decay to gravitinos |65| . the resulting gravitino abundance may give an affect on our study. 
But this can depend on the inflation model. We will also need to discuss the generation of 
/x-term and S^-term, which is related to the main decay mode of the i?-axion and other fields. 
This issue is common in gauge mediation models. (See recent gauge mediation models with the 
QCD axion [22], in which messenger fields of SUSY breaking and axion are not unified. With 

= fi'e- qA H u H d [5J and with F A ~ up to m.3/2, fJ*/B/j, problem have been solved there.) 

As future directions, we can also study another parameter region, where some dimensionless 
parameters are much smaller than of order unity. For such a case, new possibilities may open up; 
for instance, we may have much smaller masses of saxion and axino than our cases. The saxion 
dominated universe may be allowed in such cases because of the different decay properties. It 
will be also interesting to study the UV completion, where the SUSY breaking sector and the 
PQ breaking sector are unified. That will be a minimal model to solve hierarchy problem and 
strong CP problem. 
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A Supergravity effects 



In SUGRA, i?-symmetry can be broken by a constant term in the superpotential, SW = h 2 Wq. 
The scalar potential is 



^sugra — e 



where DiW = d%W + W . The gravitino mass is given by m s / 2 = /u^Wo/Mpp Here the 



K %3 D;WDjW - 3 



3 M 2 



(76) 



Kahler potential and the superpotential are given by eq.(16) and (17). Hereafter in the appendix 
we will define 

^ = X' (77) 

• i?-symmetric case: a = 1 
In this case, the vev shift is 

2A V3W x . c (rj- 2 f 2 + 2cA 2 ) 

OA — — — • — — , 0/1 — f; . 

v^Mpi Mpi 3cM 2 { 

Here Co is a coefficient of the linear term in the Kahler potential, cof 2 (A + A'), and in 

A 



order to make our perturbation valid, we need j^— <C A for /< A. The constant term 
Wq is fixed by the condition ^sugra = 0: 

W ^ -i=Mpi. (78) 
v3 

Then we find AFx = —^ 2 = —VSth^^Mpi. The axino obtains a mass 



/ 2cA 2 \ 

md = m 3/2 ( 1 + ■ (79) 

The first m.3/2 contribution vanishes when we replace the axion kinetic term, (A + A^) 2 /2, 
with ex^A+A^). This is because we can define A new = A+\og& c , where $ c is a conformal 
compensator. Hence loop effects, such that c^A + A^) • exp(^4 + A*), are important. The 
axion field A gets an F-component vev at the leading order of 1/Mp\: 

F A = c m 3/2 . (80) 

Spontaneous i?-symmetry breaking case: a = — 1 and xo ^ 
There are similar vev shifts and F A to the previous case: 

XY A c q X ( f V , A c (/V 2 + cA 2 ) (qxv~ 2 + c c) x A 

" x = ~r 1 — — 77 — x , oA - 



v^A/pi 35 \r]MpiJ u ' 3cM|j ^35 Af H ' 
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F A = c m 3 / 2 . (81) 

Here we used Wq = Mp\/y/3. Recall that we will have Fa oc Xq term, which is not 
suppressed by Mp\. The i?-axion gets a mass: 

m ^ = ^ = 2 ^/ 2 I^a J (82) 

at the leading order in xo and Planck suppressed expansions. Now the axino mass is given 
by 

™a = ™3/2 H J2 = m 3/2 I 1 + V 3r? C j2 ) " ( 83 ) 

The first m 3/ / 2 contribution vanishes again when we replace the axion kinetic term, (^4 + 
At) 2 /2, with exp(^ + At). 



B Solving mixing between X and A 

B.l Kinetic mixing in the standard form by the superfield description: axion 
mixing and fermion mixing 

Let us consider the following Kahler potential, which preserves [7(1)pq symmetry 

f 2 

K = X^X+ j —{A + A^) 2 + K X a^X{A + A^) + Kax^A + A^)X\ (84) 

As we want to solve kinetic mixing, we will focus on the vacuum in which K XA \ has the vev: 

f 2 

K -»• K = X^X + J — {A + A^) 2 + {A + A^){kX + k*X*). (85) 

Here we defined k = (K XA i)- After solving mixing between A and X, the above Kahler potential 
becomes 

/ \k\ 2 \ f 2 
K = ( 1 — X'X + J — (A + i 1 ") 2 + (holomorphic term + c.c). (86) 

A diagonalized axion A is given by 



> , kX (K x m)X 

A = A + j ¥ = A+ x y 2 T/ . (87) 



Thus when one considers the canonical normalization we find 



, A y/2A V2X , . 

k-A ^ V-V. (88) 

Ja Jy 
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Here 



With a gauge kinetic term S = ^2 + k-^j, we obtain 



Here we used K xx t ~ 1, K Aj4t ^ / 2 and K X xt-^AAt > for / > K XAt- 

Thus, with regard to / y , we can get 

J " " A; V V(^XAt)y _ fc V ™~47r M 1/2 " 
Here we used (Fx) — — \/3m 3 / 2 Mpi, (Fa) ~ ^My 2 - Thus we can find 

/a /F A 



(89) 



Furthermore, in the standard form, Fx and F-component of A is given by 



Fx ^ -(d x W)\ (91) 
Fa = -K^cW^-^^Fx. (92) 



(93) 



F .v =-tan6, (94) 

where 9 is a goldstino angle. 

Using this goldstino angle, we obtain goldstino ip x and axino a 

tx^f-fflW a^f + S f^\ (95) 
Jy Jy 

where we assumed G <C 1. In the above equation, we normalized fermions canonically and 
means the original field before solving the mixing. Thus we find (<5susyo) ~ (F|) = 0. 

B.2 Mixing between Xr and Ar 

Remaining fields to be solved are Xr and Ar. Now there is also mass mixing in the scalar 
potential in addition to the kinetic mixing: 

V = \m\ R Xl+ l -rn 2 AR A 2 R + rn 2 x A X R AR, (96) 
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Here recall that we already solved the kinetic mixing with orig inal field (X { ° ] , A$) : X R = 
X R °\ Ar = A^ + j?-Xft\ By an unitary rotation, we will obtain diagonalized fields (Xr, Ar) 
with canonical and diagonal kinetic term: 



1 2 /f \2 , L2 ( A \2 



Here 



X 



V « ^x R {X R y + -m\ R (A R y. (97) 



(o)\ 



fa 




foi\m XR -m AR \<m XA (98) 



(99) 



Here we will always have mx R > m A R because of A c /A; hereafter we will not consider the 
degenerate case. 



For Model 0, we will find 



m XA A c . ( M 2 fa ^ f, 



m\ R /M P1 V / ) fy ' fy 

Here f y ~ Mp\. 



< J -f. (100) 



B.2.1 Mixing between Xr and Ar in Model 1 



In this case, we have 

m\ R = 2^, m \ R = 2 V 2 c^, (101) 

2 2 2(-J+e)/i 4 x /a 2 / 1no x 

m xi = ^ ]^ mAii ' ^ 



.2 fa^a L , fS M l/2 ^ 

V6m 3/2 Mpi I m 5 ra^ fl 



fc^L / 1 + J 2 I ( 103 ) 



~ 2 S x 10 12 CcY 2 ( mx « Y( m& )( fa ) p 0MeV ^ 
- A8XlU GeV WcevJ UxlOSGevJ Vl0 10 GevJ V ™3/2 )' 

Here we supposed (— d + e) ~ c. Note that mass mixing in the scalar potential m XA is of O(xq) 
while kinetic mixing fy 1 is of O(xq). 

Thus 

™™ ~ k J am ' a ~ 2.8 x 10- 2 . (104) 
m x fl V6m 3/ 2M P1 
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Here we used used mx R > m A R and took f a = 10 10 GeV, m„ = 5 x 10 5 GeV and m 3 / 2 = 30 
MeV. 

B.2.2 Mixing between X R and A R in Model 2 



Recall that eq.(41): 



m XA = e K^7^ (-Qrj 2 c+ ^ ) (105) 



2// 2 \ ' A 2 , 

Cm\ (106) 



4vr M 1/2 f a 



* fc fT' aU M Cm *- (107) 

Here C = 0(1), m x = /i 4 /A 2 , and we used F A ~ m„ ~ |^ Then the magnitude of the 
mass mixing is the same order as that of the kinetic mixing: 

m 2 XA i m afa 



2 

m 



_ k aja c ^ 1Q8 ^ 

x V6rn 3 / 2 Mpi 

~ ~ 10~ 3 x C. (109) 

Here we took 777.3/2 = 30 MeV and ma ~ 10 4 GeV. 
B.2.3 Mixing between X R and A R in Model 3 

In this case, we have the similar situation to the Model 1, that is, the effect of the mass mixing 
is larger than that of the kinetic mixing: 

m 2 XA ~ " 2 ^ 2 ^ 4 (HO) 
~ - 2k J^T m x for ( U1 ) 



Then 



2fc ™ a/a (112) 

m£ V6m 3/2 Mpi 

~ -5.3 x 10~ 2 . (113) 



Here we took rria = 10 5 GeV and 771,3/2 = 30 MeV. 
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C Derivatives of F x and F A 

We need computations of djF 3 (I, J = X, X\ A, A') to obtain interactions between gaugino 
and those fields for instance. Now we will focus on the gauge kinetic term: 

l - J d6 2 SW a W a + c.c, S = ^ + k^. (114) 

As we have a gaugino mass from F A , we can read a coupling of X to gaugino pair from the 
gaugino mass at the leading order of fluctuation SfX: 

\M 1/2 = l -F A d A \og(^(S)(A)) (115) 
- k^^jf(d x W)^ (116) 
= \(M 1/2 )(l + 5fX(d x \og(K XA ,))) + ... (117) 



^(M 1/2 )(l + ^) + .... (118) 



Here we assumed K XA \ is a polynomial in X, so we can gain (d x log (K XA ^)) ~ V(^0 un l ess 
we have a small X vev or cancellations like Model 2 below. We used a notation (Mi/ 2 ) to 
distinguish dynamical fields from a parameter. 

C.l Derivatives of F A in Model 

~ -d x ,F x ~ -m 3/2 , 0^ = 0^*-^-^, (H9) 

^ ~ ^(^Va, O^-^A, d A F A = d A <F A ~m s/2 . (120) 
Mpi V m x / Mpi 

C.2 Derivatives of F x and F A in Model 1 

d x F x = d x <F x ~F A xl^+(£j 2 y d A F x = d A ,F x ~F A x (£) 2 . (121) 
d x F A = - — , d x ,F A ~F A x , d A F A = d A ,F A = ma (122) 

Xq 
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C.3 Derivatives of F x and F A in Model 2 



6 X F* ~ ^(!)V(£f< ^ = -e^3^--V2A^i (124) 

9 A F A = 8 Ai F A = -C 3 F A ~ -m g . (125) 

C.4 Derivatives of F x and F A in Model 3 

9^ = ^-W^Y, S xt F* = -^~WZY (126) 



A \eAJ ' A V A , 

d^F x = S^F* = k2 j§^ pA ( 127 ) 
SxF 4 = ^xti^ = ^ J ^FA = 5 At ^ = m a . (128) 

D Some decay modes 

Here we will exhibit several decay modes. The results can be rough and just show their order 
of magnitude. 

D.l The Polony i and the i?-axion decay 

• T(X R -> 2Xi) in Model 1 

Through an interaction A Xj Xi in the Kahler potential K = A 2 X^ X, we can 
obtain 

nXn^X.X.) * (129) 

. T(X -> 2y 3/2 ) 

With an interaction K = X 2 or K = A 2 X^X for the 72-axion in Model I 19 we can 

find 

I TYv 



We have an interaction between goldstino ipx and i?-axion X/ through K = A 2 X t X: ^J^d^Xiipxa^ipx. 
A contribution from KA~ — - — j-* 1S a w iU vanish because of the vev (X) ~ xq. 
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F(X -> V3/2 + o) 

Through an interaction Kax, we obtain 



r( ,^, 3/2 + a) K (f )^^. (13!) 



rpr -> 2 ) 

Through an interaction with mixing between A and X, 



£ f + ) , (132) 

OVT V f a Ja Jy fa 



we obtain 



For Model 0, we will not have a loop factor. 
F(X -»• 2a) 

Through an interaction 

X 

nia—j= ad for Model 1,3, (135) 

ma^-dd { 1 + — f ^ ) I for Model 2, (136) 



kfy \ 2 V m, 
we obtain 

r(X -»• 55) ~ ( ^ ) for Model 1,3, (137) 

2 



) 1 + T(^) M°del 2,0. (138) 



Note that (X) is a dimensionful vev of X and the above interactions are obtained via 
V 2 A 2 I{A + Ai) 2 X*X and + At) 3 . 

F(X -»• 2A) 

For a mode X — >■ AA, dxF A or d x ^F A becomes a coupling of X to gaugino pair: 

N„ fM l/2 \ 2 

T(X^AA) ~ -J-mxi^yj for Model 1,3, (139) 

jrnx — ^ for Model 2,0. (140) 



167Tfc2 V /, 
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M 



Thus a partial decay width of this mode can be suppressed by ( "*2.{ 2 ) for Model 1 and 



2 



-2 



3 or (l|r) ( 1 + ^ (~^f)~ ) for Model 2 and °> comparing with T(X -»■ 55). If we 
have if = c^(A + A^)H u H d + cb(A + A^) 2 H u H d , T(X — > h u h d ) would be comparable to 
this mode, supposing ~ 1CT 2 , eg < 10~ 4 and \i ~ M^- 



r X/ in Model 1 

fi n rl + Vi twit rrVi qti i nf or a r>\- 1 r^n — i . 



For Model 1, we find through an interaction —i M i /2 XX + c.c. 



N n ( Ml /O 

r(X/ AA) = ^m Xl ' 
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Ax J \ 



4M l/2 

1--^- (141) 



For < 2Mu2, the channel X/ — > XX is closed, then Xi — > bb can become the main 
decay mode: 

r(X/ ^ 66) " x A ml-mlj f ~ mf/ (M2) 

For a mode of Xj —> tt, mj, should be replaced with m T . 
T(X R -> 2$) 

We denoted MSSM scalar fields as For decay mode X R -> 2$, K = k (^) 2 (A + 

2 2 

7 4t) 2 <j)t<|) j s important. These decay occur via soft scalar mass ^jXr&$> or ^Xr&Q. 
These amplitude is given by 

T(X R -4 2$) ~ £ ( Jp-) 2 (£SL) 2 x T(X - AA). (143) 

iVg VMl/2/ WW 

Here JV m is the number of these decay channels. Supposing we have K = c^(A + 
A^)H u H d + cb(A + A^) 2 H u H d and ~ m 2 ,, T(X -> H u H d ) from this interaction 
would be comparable to this mode. The decay may also occur via derivative interaction 
k (^) 2 (A ) jf&d 2 $ [66J. (Note that A includes (A ).) However the effect of this 

interaction is much smaller than the above result. For Model 0, we will have similar result 
to the Model 2. 

T(X R -»• 2a) etc. 
Through an interaction 
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we will obtain 

1 111 Till 

T(X R ^aa) * j^^. (145) 

Note that the above interactions are obtained via rj 2 K 2 ^(A + A^) 2 X^ X and ^r L (A + A^) 3 
for Model 1, 3 or via ^-^-{A + A) 3 for Model 2. For mx S> mi R , we have the same 
partial decay width of Xr — >■ A^A/j as this. 

D.2 Saxion decay 

T(A R -> 2a) 

Through an interaction Cs^d^ad^a in the Kahler potential g |f (vl + A^) 3 , we can 
compute 

r <^ + = sSS^lf 1 - (146) 

r(^ R -> a + X/) 

From the previous interaction, we have 

r(A fl -> v 3 /2 + a) 

Through an interaction K = r/ 2 | (A + A^-X^X, we can find 

r ^-^3 /2 a) - (148) 

r(A fl -> 2V> 3 / 2 ) 

Through a mixing between Xr and .Ar, we can find 

r « i ^*'*'-» = (149) 

But this will be suppressed by (y 2 ) , compared to T(Ar — > V3/2 + a)- 

r(A R -> 2 5 ) 

Through an interaction 

^jrW" ( 15 °) 
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we obtain 



N n ( a \ 2 m 6 A 

r(A^ S9 ) = -f- (151) 



For Model 0, we do not have a loop factor in terms of /. 

T{A R -> 25) 

Through an interaction 



C^nia^aa (152) 



in ^^(^4 + A" 1 ") 3 , we will obtain 




n^R^aa) ~ -« Jl-4-^-. (153) 



2 



Here we used 8a r F a ~ 2ma. Note that we may have also an interaction rj 2 jk m ajf aa 
via KA~ 2 = —i] 2 ^(A + A^) 2 X^X, where d' is a coefficient of mixing matrix Uu, e.g., 
<f = T] 2 (d — e) for Model 1. But this is at most comparable to the above interaction. 

F(A R -> 2A) 

Note a decay mode Ar — >■ AA comes from a gaugino mass interaction ^m-a^AA with 
using 8 Ar F a ~ 2m a . Thus 

m * (154) 




2 

mi. 



is suppressed by N g (^) , compared to r( J 4/j — >■ da). Note that r(^R — > h u hd) would be 
comparable to this mode, supposing we have K = c^{A + A^)H u Hd + cb(A + A^) 2 H u Hd, 
c ii ~ i^- For Model 0, fc^ma should be replaced by m 3 / 2 ~ 

r(A fl -> 2$) 

For a decay mode — >■ 2$, K = k (|^) 2 (A + Af) 2 <J>t<£ is important. This decay can 
occur via soft scalar mass j^ m J^ Ar&$>. The amplitude is given by 

F(A R -> 2$) ~ ^ f-^V x r(^ -> AA). (155) 
iv s \mA R J 

The decay may also occur via derivative interaction (|^) 2 (A)) j^&d 2 ^. However, this 
effect is much smaller than the above result. For Model 0, we will have a similar result, 
except that an interaction ^r/iQ$t$ becomes relevant. 
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T{A R 2X T ) in Model 1 

Note that we have -j^- ~ rj 2 for Model 1. (Recall that (X) is a dimensionful parameter.) 

d(X) 

Then we have interaction (d fl X I d^X I + 5d^,X R d^X R ) for Model 1. Thus, through 

this interaction we can find 

riAR^XjXj) ~ g^(^) for Model 1X 7 . (156) 

Here ^ ~ rj 2 for Model 1. 
D.3 Axino decay 

Through an interaction in gauge kinetic term 

-^-aa^X a F^ + c.c, (157) 

the heavy axion mainly decays into gluino and gluon because of the strong interaction. 
Then we have 

F(^A +9 ) * iV^I (158) 
For Model 0, there is no loop factor in terms of /. 

r(s^^ + $) 

Here we denoted ip as the MSSM matter fermions. Note that we also have the fermion- 



sfermion- axino interaction 



° 2 ^logf^l^a (159) 



V2vr 2 fa b \M lf 

via loop correction by the MSSM gauge interactions [67]. However, this interaction is 
irrelevant since axino mass is large: amplitude T(a — > ip + $) will be suppressed by 
^ (a log (i7^)) (M 1/2 /m 5 ) 2 compared to r(5 -»• A + 5 ). 

r(o ->• V3/2 + «) 

Through an interaction 

aa^^^dua (160) 



Mpi 
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or equivalently, 




this decay will occur. Here the above goldstino interaction originate from — 7] 2 A 2 ^(A + 
A^) 2 \X\ 2 or ^(A + At) 3 . Then we find 



T(a ^3/2 + a) 




(162) 



• r(a ^3/ 2 + X/) in Model 1 

In Model 1, we can have an i?- breaking interaction ma^j^^aipx in ^AX- Thus we find 
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